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NAG8-684: Process modelling for SpaceStaUon.rperlmeot, 

Annual Report for the period 10-1-87 through 9-30-88 

performed during the first year 10/1/87-9/30/88 
involved analyses of crystal growth from the melt and from solution 
Thf particular melt growth technique under investigation^ 
directional solidification by the Bndgman-Stockbarger me _ 
types of solution growth systems are also being studied. One 
Xth from solution in a closed container, the other concerns growth 
S protein crystals by the hanging drop method. FoUowing discusmons 
witt. Dr R. J. Naumann of the Low Gravity Science Division at MSFC 
decided to tackle the analysis of crystal growthf^m the me, 
earlier than originally proposed (see Task Schedule, section 2.3. pag^ 
12 of Proposal). Rapid progress has been made in th s • 
on schedule and full calculations have been underway for some • 
Progress has also been made in the formulation of the two solut o 
growth models. More details are available in the Mowing sections. 

Numerical methods and code development 

For all three model systems PHOENICS. a commercially avatlable 
finite volume code, has been used extensively. In addition, a code as 
been developed, which employs a Chebyschev-collocation pseudo- 
spectral method, and applied to both steady and time-dependent cases 
for the analysis of crystal growth by the Bridgman-Stockbarger 
technique. This method has been used to verify our ear l er 
calculations, increase accuracy (particularly in the vicinity of the sohd- 
liquid interface) and to extend the range of parameter space under 
investigation. We are currently using this code, and are also processi g 
L results of calculations carried out over the last few months^ 
Preliminary results for thermosolutal convection have also bee 
obtained and will be extended. Further calculations are also planned 
which will cover a wider range of parameter space and involve more 

realistic (g-jitter) disturbances. 



2 


Crystal growth from the melt by the Bridgman-Stockbarger method 

The formulation of the model was completed during the first 
semi-annual period and progress in this area is advancing faster than 
originally anticipated. Preliminary results were reported in detail in 
the first semi-annual report and a publication is currently in 
preparation describing the first year’s results. 

The goal is to examine the effects of the microgravity 
environment on the compositional uniformity of crystals grown by the 
Bridgman-Stockbarger technique. In past space experiments this 
technique has been used for the directional solidification of crystals, 
and it is anticipated that this technique will be used in the Materials 
Technology Module (MTL) and the Materials Processing Facility 
(MMPF). During the first year our modelling efforts in this area have 
resulted in the development of two- and three- dimensional models of 
directional solidification which have been used to determine the 
sensitivity of crystal compositional uniformity to both time-dependent 
and steady residual accelerations. Selected results are described in 
detail in appendix A. 

To date we have shown that: 

a) A steady background level on the order of 10‘® - 10'5 g can be 
tolerated provided that the acceleration vector is aligned with the 
axis of the growth ampoule, and provided that no accelerations with 
frequencies less than 10' 2 Hz (and amplitudes on the order of the 
steady component) are present. 

b) For a fixed growth rate, the amount of lateral segregation is 
very sensitive to the orientation of the steady component of the 
residual gravity vector. The worst case appears to be when the 
acceleration vector is parallel to the crystal interface. At growth rates 
on the order of microns per second, this orientation leads to 
compositional non-uniformities of 22% when the magnitude of the 
acceleration is V 2(10)'6 g. If, however, the growth rate is lowered by 
an order of magnitude, the non -uniformity is reduced significantly 
(down to 4-5% in this case). 
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. response of the solute field, and the lateral non-uniformity. 

frequencies above 1 hz , f 10 -6 a in addition, 

4 . i n -3 Hy at amplitudes on the order ot iu g. 
response at 10 Hz P 4 ,r_„ c r»f a steady component 

additive effects were Le to 

and a low frequency comp f jti j n concentration, 

significant lateral and longitudinal non-unifonnlUes 

« r it r— , » - 

" Ise For exate a pulse with a one second duraUon. or a 
the impulse. For examp P ^ effect 0 n the segregation 

combination of snchpuse ^ ^ ^ nature of the response 

levels at pulse amp u d durat ion of the impulse, 

depends on the magnitude. direcUon an A s0 . c alled 

and whether sequential opposing impu s comp letely offsetting 

■compensating" double pulse will result ing 

effects. For the case we examined howe ^ ^ sequentlal 

cortm^ thTl" — Further investigation of more 
^tieCresTnecessary £££ 

r::z ^ to >». ,mpo rt ant 

consequences for transient behavior in crystal growth systems. 

It should be borne in mind that to date - — ons havemfiy 

C Td ed th a at P tr1 gill g ^e P vT” e re r sX a aTacceleraUon the amount of 
Tamral segregation can ^^rel^ 

Schmidt number lequl^lent to inc distributlon 

decreasing solute diiiusivuyi 
research program. 
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Thermal convection In a square cavity 

This work was carried out in order to investigate the reliability of 

order of magnitude estimates of the effects of residual accelerations on 

transport conditions in simple fluid systems. 

The chosen systems involved thermal convection in square and 
rectangular cavities. The rigid upper and lower boundaries were 
maintained at constant but different temperatures. The temperature 
differences examined ranged from 5 K to 50 K over distances of 1-5 
cm. The sidewalls were insulating. The applied residual gravity vector 
was time-dependent in both magnitude and orientation. Time 
dependent and steady accelerations were examined. It became readi y 
apparent that while the fluid responded to the body force. e 
temperature field was relatively insensitive. At Prandtl numbers ° 
order one the temperature field exhibited minor oscillations about the 
basic diffusion profile (a simple linear profile for this case). However 
none of the cases examined showed a significant fluctuation in heat 
flux at the low temperature boundary. The major reason for this 
appears to be the strong influence of our thermal boundary conditions 

on the heat transfer rates in the system. 

In order to check the reliability of order of magnitude esUmates 

the g-sensitivity predicted by the Langbeln-Tlby approach was also 
tested. The basic scenario is that described in their paper HI. In 
many fluid physics experiments, when particle motion is observed or 
the behavior of a growth front is studied, one can assume that dr/r [the 
average density change] in the test volume has the order of magnitude 
10-2 -the diameter of the test volume is chosen to be 2 cm, 
temperature differences of 5 K are applied and an accuracy of ST 
0 05 K is required ..." In addition, we assumed adiabatic sidewalls. 

The estimated tolerable g-level for the the above requirement is 

SiVel aLax = 2[ 0)2 + (v/ L 2 I 2 ] 1 / 2 [ oi 2 + (x/ L 2 l 2 )!/ 2 . 
and is based on a "typical velocity" magnitude given by 
O l u(r)] = (8p/p) a 0 co _1 l co 2 + (v/ L ) 1' 
where co is the frequency of the residual acceleration, the kinematic 

viscosity of the fluid is v = 10" 2 cm 2 s’ 1 , k = 10 cm s is t e 
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thermal diffusivity, L is the characteristic length (taken to be 2 cm) 
and a 0 is the magnitude of the acceleration. The predicted order of 
magnitude of the velocity u for a steady acceleration of 10 -6 g is three 
orders of magnitude higher than the maximum velocity computed 
from the solution of the full non-linear equations. In every case we 
examined, the maximum velocity estimates were at least two to three 
orders of magnitude higher than our maximum computed velocities. 
As a consequence, the effect of residual accelerations tend to be 
overestimated. These discrepancies can be explained simply in terms 
of lack of a suitable length scale on which to base the estimates and 
the fact that order of magnitude estimates cannot, in general, properly 
account for the geometrical and thermal nature of the boundary 
conditions. 

Effects on transport conditions during protein crystal growth via 
the hang ing drop method 

The s tar ting point for this work is the simplified one-dimensional 
diffusive model of Fehribach and Rosenberger [2]. To examine the 
effects of residual acceleration on this type of crystal growth technique 
buoyancy driven convection must be included in the model. In 
addition, the nature of the interfacial conditions between the hanging 
drop, the "wall" to which it is partly attached and the vapor with 
which it is in contact, result in a complicated moving boundary 
problem which cannot be adequately handled using our available codes. 
Hence, development of a suitable numerical algorithm is in progress. 
The method under consideration involves a front fixing method 13], a 
finite difference approximation in space and an alternating direction 
implicit (ADI) method. 

Transport conditions during crystal growth from solution 

Triglycine sulphate crystals were grown from aqueous solutions on 
Spacelab 3 using the Fluid Experiment System (FES). Under the 
assumption of zero-gravity, Yoo et al. [4] have modelled the diffusive 
transport conditions for this system and estimate the influence of 
convection using order of magnitude considerations. Their results 
suggest that convective effects cannot be ignored. 
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in order to obtain more quantitative information concerning the 
effect convection we have formulated a model based on that of Yoo 
et al 141 but have included momentum as well as mass a 
Insfer equations. As a practical measure the increase in computation 
consequent to the introduction of the momentum equabon has bee 
offTsornewhat by simplifying the geometry of the computational 

domain. 

C ° m c<mnmnicaUons with MSFC personnel (Dr. R. J. Naumann and Mr. 

F Rutledge) have taken place on a regular basis. A detailed semi-annual 
report was submitted in April 1988. Specific information wa 
supplied to Dr. R. J. Naumann (MSFC) concerning micrograv ty 

requirements to be used in testimony to Congress. 

Selected results of our analysis of the Brldgman-Stockbarger 

technique have been the subject of severa. presentations mcludmg the 

MSFC Distinguished Lecture Series, May 31*t 1988 (given y 

SSL. «o««. 

Science and Space Experiments session at th 

Conference on Drops and Bubbles in Monterey California Septemb 
19-22 1988. 

Planned Work for the third six month period 

During the first six months of the second year the analysis of 
crystal growth from the melt will be brought to completion. Wor 
remahiing to be done in this area includes a more detailed 
examination of thermosolutal convection and an examination of more 
"ic residual accelerations (g-jitter). During this period the 
preliminary stages of models for crystal growth from ^ution (c os 
container) will be completed and full scale calculations will begin. 
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APPENDIX A 


Response of convective-diffusive transport to spatial 
and temporal variations in effective gravity 


By J. IWAN D. ALEXANDER JALIL OUAZZANI 
AND FRANZ ROSENBERGER 
Center for Microgravity and Materials Research 
The University of Alabama in Huntsville, Huntsville Alabama 35899 


The reduced gravity environment on board a J P nhvsl^sts the 
earth orbit gives materials scientists and ^ JL at reduces 

nnnortunitv to undertake experiments under conditions that reduces 
or Eliminates buoyancy driven fluid convection in comparison toearUi 
Cased Auditions As a consequence, the relative importance of heat 
and mass transDort by diffusion is increased. In crystal growth this 
mitfht be exoected to lead to a more uniform crystal composition than 
™Sd be obtained under terrestrial conditions. The process of^rysfcd 
growth by the Bridgman technique is chosen as a case study . TVo 
toree dimensional numerical models are used to examine the 
resDonse of heat, mass and momentum transport to condition 
characteristic of the microgravity environment. It is shown that the 
orientation of the experiment with respect to the steady component of 
the residual gravity is a crucial factor in determining the suitability o 
the spacecraft as a means to suppress or eliminate unwanted effects 
caused by buoyant fluid motlom The process is also extremely 
sensitive to transient disturbances. For example, a 3x10*3 g impulse of 
one second duration acting parallel to the interface of a growing 
cwstll produces a response to the solute field which lasts for nearly 
2000 seconds. Consequently lateral and longitudinal compositiona 
variations occur over a length of nearly 6 mm in the grown crystal. 
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1. Introduction 

For over two decades, particularly since the operation of Skylab in 
the early seventies, the prospect of undertaking experiments aboard a 
spacecraft in low earth orbit has appealed to scientists from several 
disciplines for a variety of reasons. In addition to serving as useful 
base for cosmologists and astronomers, an orbiting spacecraft 
provides an opportunity to study physical phenomena in the absence of 
the 9.8 ms- 2 (1 g) acceleration experienced in an earth based 
laboratory. The various phenomena that have been the subject of space 
experiments have been described in a number of recent books and 
articles (Abduyevsky 1985; Feuerbacher et al. 1986; Hazelrigg & 
Reynolds 1986; Walter 1987). 

Physical systems which respond to buoyancy forces will exhibit 
different behavior in space than in an earth based laboratory. A simple 
example is the behavior of bubbles or liquid drops with densities that 
differ from the host liquid. At the earth's surface, these move up or 
down in response to gravitational acceleration In a spacecraft they 
will no longer exhibit such predictable behavior. More complex 
examples include the surface shapes of liquid films, bridges and drops, 
and mass transfer in crystal growth systems (which we shall examine 
in this work). Under earth based conditions mass transfer conditions 
can be strongly affected by buoyancy driven convection in the (fluid) 
nutrient phase. In the reduced gravity environment of a space 
laboratory, mass transfer is more likely to be controlled by diffusion 
(Rosenberger 1979; Sekerka & Coriell 1979, Hurle et al. 1987). 

As the number and level of sophistication of space experiments 
has increased, our understanding of the characteristics of the space 
laboratory environment has also changed. What was initially referred 
to as "zero-g" has become "micro-g" and more recently transient 
acceleration measurements revealed "milli-g" perturbations caused by 
crew activities, machinery vibrations, attitude changes etc. (Chassay & 
Schwaniger 1987) Indeed, it has become clear from measurements of 
the acceleration environment in the Spacelab (Chassay & Schwaniger 
1987; Hamacher et al. 1987), that the residual gravity levels on board 
a spacecraft in low earth orbit can be significant and should be of 
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concern to experimenters who wish to take advantage of the low 
gravity conditions on future Spacelab missions an o 
nlanned Space Station. While accelerations may be orders o 
magnitude lower than that experienced at the earth's surface, they are 
nonetheless finite and. thus, pose potential problems for certain types 

0f To 6 d^r^alyses of the effects of a low gravity entdronment have, 
with a few exceptions (Kamotanl et al. 1981: Chang Brown 1983 
Polezhaev 1984; McFadden et al. 1988) been restricted to either 
idealized assumptions or order of magnitude estimates (Spradley et al. 
1975- Camel & Favier 1984. 1986: Langbein 1984. 1987; Boudreault 
1984- Monti & NapoUtano 1984; Monti et al. 1987). The validity o e 
Zus estimates has only been demonstrated for a few special cases 

tR0U ^ U a case stady. we examine the effects of time-dependent and 
spatiad variations ln the effective gravity vector on solute — ution 
in a simplified model crystal growth by the Brtdgman-Stockbarg 
technique. In section 2 we discuss the phenomenon of segregate 
which when influenced by unsteady and spatially varying transport 
conditions ln the nutrient, leads to undesirable compositional 
variations. In section 3 we discuss the basic heat, mass an 
momentum transport model of the Brtdgman-Stockbarger technique 
•Hie results obtained from our numerical models of the effects 
residual steady and time-dependent acceleration during the melt 
growth of a doped semi-conductor crystal are presented in section 4 
In d discussed in section 5. Our results show that for this sped 
directional solidification system, residual g-levels of g are 

adequate only If an optimal alignment of the system with respect 
the residual gravity vector is maintained. 


1 2 Solute Redistribution 

Consider a two component solution system. The equilibrium 
between the solid and liquid phases ln such a system can e 
represented by a binary phase diagram (Chalmers 1977). xcep or 
the case of either a pure or a congruently melting material, the 
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tol Uquid with which It is in equilibrium. Such compositional changes 
on solidification occur in the majority of systems of interest to th 
mater!^ science community, particularly those materials needed for 
device fabrication (Rosenberger 1979). During ^ £ ay*£ £» 
melts containing two or more components this ^ 
rise to a phenomenon termed segregation or redistribution. 

Segregation is of great importance for many aspectsofmaterial 

dT rforexampl! 7 me" puXtion of Tatertals or the 
SUS composition of the — ^ " 

!;,ow the liquid phase to adjust its composition throughout in 

accordance with the equilibrium phase diagram. This is ee * use 
accordance 4 Uquld are typically slower than the 

rr.ficaUon rate. As a consequence. 

persist at the phase boundary, compositional gradie 

established in both the solid and Uquld phases (see Fig. 2). For case 

which local equilibrium conditions prevail at the ph^e bou^ 
the interface! transfer of a component may be described _ta tenns^of 
equilibrium distribution coefficient ^-nberger ^ ^ 

gfven “cT/Tm\ whe^e e s and cm are respectively the solute (B) 
concentrations in the solid and liquid. The value of k reflects whettier 
the solute is preferentially rejected or Incorporated upo 
solidification. The equilibrium distribution coefficient depends orfiy o 
the thermodynamic equilibrium properties of the system. 

independent of the mass transfer kinetics. 

in order to predict the operating conditions under which 
desirable dopant uniformities can be obtained, for example for doped 
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semi-conductor crystals, it is 

between tire - which the 

:irCrcUi -fleets the local J- — 

rn“rmCn ^ ^ ~ «“£ ^ 

result in compositional non-uniformities in convect i 0 n in the 

M been recognized for ' som^ ^ composltlonal 

nutrient phase can, in some * crystals (Pimpuktar & 

variations in melt grown semi-con u these 

Ostrach 19811 MOller 1982: MMr - 

cases, for example growth by the Bong attractive 

- — e 'low^l^^e^^avliy^e^dtonment afforded by a 
feature of the low eitecti g ' reducing or eliminating 

spacecraft, is the possibility o su U ng diffusion controlled 

buoyancy driven convection, mid ^ L favor compositional 

conditions (Hurle e • . Uowlng mo del of crystal growth by 

r^Tal^gZ^Ser technique ^—te.e 

under steady ume ' 

dependent reduced gravity conditions. 


rild'lg adiabauc sidewalls. “^"01^ 

— r;: modced by 

“ pS Xed me y it Of dilute 

velocity V„ at the top of the computational space both 

- L'rrrsr.?, — - . — 
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at a distance L from the top of the computational space. The 
temperature at the interface is taken to be Tm. the melting 
temperature of the crystal, while the upper boundary is held at a 
higher temperature Th. In an actual experiment, owing to the finite 
length of the ampoule there is a gradual decrease in length of the melt 
zone during growth. In this model transient effects related to this 
change are ignored. Thus, it is assumed that the ampoule is 
sufficiently long for these effects to be negligible. The only transient 
effects to be considered will arise directly from the time-dependent 
nature of the residual gravity field. We also assume that the 
contribution of the solute (dopant) to convection is negligible. 
Convection is driven only by thermal gradients. It has been shown that 
curvature of the solid-liquid interface can result in significant lateral 
compositional non-uniformity (Coriell & Sekerka 1979; Coriell et al 
1981). We wish to focus attention on the influence of convection on 
the composition of the crystal and thus, have chosen to examine the 
planar interface case. For the two-dimensional model the dilute binary 
melt is assumed to occupy a rectangular region £2 which is bounded by 
planar surfaces (see Figure 3). In the three-dimensional model the 
rectangular region is replaced by a circular cylinder with a diameter 

equal to L. 

The governing equations are cast in dimensionless form using L, 
the effective melt length, K/L (where K is the melts thermal 
diffusivity), PmK 2 /l2 , Th-Tm. and Coo to scale the lengths, velocity, 
pressure, temperature and solute concentration, respectively. The 
dimensionless equations governing momentum, heat and solute 
transfer in the melt are 

+ (grad u)u = -gradp + PrAu +RaPr0g(t), (1) 
9t 

div U = 0, 


— + u. grade = A 0 , 
dt 

— (— + u* gradC) = A C, 
Pr 3t 


(3) 


( 4 ) 


7 


where u(. t) represents the velocity. 9 = (T(*.t) - T M ))/(Th - Tm) the 

the hot zone and the crystal Interface) and C tePresenU 
concentration. The parameters Pr=v/K. Ra= gP(T H -T M )L /KV ana 
S c=v/D are respectively the Prandtl, Rayleigh and Schmidt number^ 
The term g(t) In equaUon (1) represents the (time-dependent) gr ty 
^ctoT^he value of g m Ra Is taken to be 980 cm s- l.e. equal to ffie 
terrestrial acceleration. Thus, the magnitude of g represents the ratio 
between the actual residual acceleration and g. Table 1 lists the orms 

n f j»rt) W e have examined to date. , _ 

The following boundary conditions apply at the crys 

interface ^ 

9 = 0 , 

( 6 ) 

u-N = Pe/o, 


N x u x N = 0, 


ac 

dz 


PeSc 

Pr 


(1-k) C, 


(7) 

( 8 ) 


tzzsxs srx-.xsxxss:. ' 

= p M /p 9 and k is the distribution coefficient. We define the measure of 
compositional non-uniformity in the crystal at the interface to e e 
lateral range in concentration given by 

(C 9 max- C smin) xl00% . 
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where c s is the (dimensional) solute concentration in the crystal, mid 
r,s L average concentration. At the "Inlet'' (z-0) the following 

boundary conditions are applied 
9C PeSc 

dz 


Pr 

0 = 1 , 

u-N = Pea, 


(C-l), 


(9) 


( 10 ) 

(ID 
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( 12 ) 

N x u x N = v, . 

fo\ a fQl exoress conservation of mass at the crys me 
Equations (8) “ d Uvely . Equations (6) and (10) guarantee 

conditions are applied 


gradC 


e = 0, u-N =Pea, 

c w ’ 


e • u = 0, 


W 


along with 

0 = 1 , 

in the isothermal zone and 
grad 0 • e„, = 0, 


(13) 

(14) 

(15) 


w 


▼ ▼ 

m the adiabatic zone. Here e w is the normal ^ ^^ mace 
While the above model does not sWcOy PPV ^ walls 
(for example details of the heat trans er -generic" model with 

neglected), it nonetheless serves ® soUdiilcation 

'EZ'SZZL Of the low gravity 
en “.:~ are limited to 

corresponding to dUute gallium-doped germam less group s 

thermo-physical properties and the associated d cakulatlons 

and operating conditions are given ta " a ^^° erence T „ - t m 
the length, L, of the furnace and the temperatu 
were taken to be 1 cm and 115° C. respecUvely. 

2.2 Method of solution th code phOENICS 

The governing ^^TdU wse) PHOENICS embodies a 
(Spalding ISSH Rosten P formula0on (Patankar 1980). It 
finite volume or finite ao orevious section 

represents the governing equations the 

as a set of algebraic equations. These equations v 
consequence of integrating the differential equation over the finite 
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volume of a computational cell (and. for transient, problems over a 
finite time) and approximating the resulting volume, area and time 
averages by Interpolation. For the 2-D calculations discussed In this 
report we employed a 40x39 grid. The 3-D calculations were 
performed In a circular cylindrical domain with 20 nodes In the radial 
direction. 12 in the azimuthal direction and 39 in the axial direction. 
We found that for the 2-D calculation less than 40 points in the 
direction parallel to the interface resulted in poor convergence of the 

somtefleldheme emp , oyed has toe same accuracy as a finite difference 
scheme which Is of the order Ax, where Ax is the distance between 
the grid nodes. The time scheme is implicit and thus unconditionally 
stable. However, small time steps are required to obtain accura e 
solutions. For highly non-linear flows the use of under-relaxation is 
necessary to eliminate divergence and to ensure good convergence of 

We compared our results (to be discussed in detail in the 
following sections) with the results of Chang and Browns 
axlsymmetric calculations. Our 3-D axlsymmetrlc computations were 
found to be in agreement with their work. In addition we carried out 
full 3-D (non-axisymmetrlc) steady calculations in order to c 1 e 
our 2-D results. For one set of examples we found that with ^2(10)- g 
parallel to the Interface the percentage compositional non-uniformity 
predicted by the 2-D calculation ft = 152%) was 50% higher than 
predicted by the full 3-D calculation. This difference may be 
attributed to the increased surface to volume ratio which increases the 
effect of the "no-slip" boundary condition. The presence of the rlgi 
walls somewhat retards the flow and. thus, in comparison to the 2-D 
case the degree of solute redistribution is decreased. A comparison 
between 2- and 3-D calculations for -j2(10)-6 g revealed no significant 
difference in At this magnitude of the gravity vector, convection is 
localized and weak. The increase in surface to volume ratio for this 3-D 
case has little influence on the transport conditions. 


Results 
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3.1 Steady accelerations: 2-D 

The temperature field (Figure 4.) was found to be Insensitive to 
the slow convective flows. This is due to the low Prandtl number of 
the melt. Note that there are lateral as well as longitudinal 
temperature gradients in the system. The consequent density 
gradient is responsible for driving convection in the melt even when 
residual acceleration magnitudes are as low as 10- 6 times that 
experienced under terrestrial conditions. 

The amount of lateral concentration non-uniformity varied with 
both magnitude and orientation of the residual acceleration. When no 
residual acceleration is present the isoconcentrates are all parallel to 
the planar crystal-melt interface. Figures 5-8 illustrate the solute and 
velocity fields for two cases where the residual acceleration is parallel 
to the crystal-melt interface. The values of lateral non-uniformity % for 
these and other cases are listed in Table 3 along with the orientation 
and magnitude of the associated acceleration vector. For a given 
magnitude, accelerations oriented parallel to the crystal interface 
resulted in the maximum compositional non-uniformity. 

For the purpose of comparison with results of calculations 
corresponding to terrestrial gravitational conditions we refer to the 
results of Chang and Brown (1983) for growth under axisymmetric 
conditions (i.e. gravity is parallel to the ampoule axis). Their results 
show that the amount of compositional non-uniformity in the crystal 
varies non-linearly with increasing magnitude of the gravity-vector. 
The maximum non-uniformity occurs at 10' 2 g, an intermediate value 
between the value of gravitational acceleration experienced at the 
earth's surface, and the 10* 6 -10- 5 g values characteristic of the quasi- 
steady component of the effective gravity in a spacecraft. Unless the 
acceleration is lowered (with respect to terrestrial conditions) by four 
orders of magnitude, the compositional uniformity may not be 
improved. Similarly, we have evaluated the conditions under which 
lateral non-uniformities are a maximum when the residual gravity 
vector is parallel to the interface. In this case, for the particular 
system under consideration, the maximum lateral non-uniformities 
will occur for accelerations of the order 10* 4 g 
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80 parallel to the crystal-melt 
were fully three aunc . t was found to be 

interface. The compositional non-uniformity^w^ ^ ^ ^ 

approximately 10 % J^Cried out forll *11 = 10 ‘ 4 . 10‘ 3 

D analogs. These calc tal are shown in Figures 9 

and 10-2. The isoconcentrates in £ ampoule axis . The fully 
and 10 for a section cut perpen c t/ 2 ( 10 )- 6 . At the 

" — r rs sH t—*- 

higher value of th sJoflOl ' 6 g E = 26%. Thus, for 

rrrirr •xzjzzx . ... » 

predictions. 

3.3 Time-dependent acceieratlo-i frequ ency 

3 . 3.1 Single frequency; combined steaay 

A number^ d^erent types of periodic disturbances were 

examined. Single frequency disturbances of form *> - g 

cos( 2 iCC 0 „t) were examined with go ' crystal . me lt 

s-. -; e 

10-2. 10-1. 1 and 10 Hz. For frequencies ^ ^ £d. 

were no discemable effects on the solute flelds^ The vetocity 

. tn the oscillatory disturbances. For the case ui 
however, respond to the oscilla ly significant. 

Hz (at 5x10-6 g) the response of the 15% were 

Lateral and longitudinal " ^-uniformity as a 

" d of uSTand hTghlight the addiUve effect of oscillatoiy and 

Stea TOe°^fect e o^a°nmUi^l^r^uency r dUturbance is illustrated in 

i x, has three components, consisting 
Figure 13. The acceleration has three co P 

steady and periodic contributions and has the form 
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g(t) = go + tfi cos(2tc10'3 t) + g 2 cos( 2 ti 10'2 t), 
where 1 1 go 1 1 * V2(10)-6, ||gill = 3^2(10) -6 and 1 1 g 2 II = 3V2(10)- 5 . 
The magnitude of £ is seen to vary with the frequencies of the 
acceleration. 

3.3.2 Impulse-type accelerations 

Four cases of impulse-type disturbances were examined. All 
impulses were superimposed onto a steady background acceleration of 
V2(10)' 6 g which was oriented parallel to the crystal-melt interface. Of 
these the one second duration pulses had the most dramatic effects. 

Figure 14 depicts the flow field immediately after a one second 
3(10~ 3 )g impulse oriented anti-parallel to the background acceleration. 
Figures 15-20 illustrate the development of the solute field following 
the impulse. Note that the effects are long lasting. The velocity field 
relaxes back to the initial state after some 300 seconds. The response 
of the solute field lags behind. The effect of the impulse is to initially 
re-orient the flow field (compare Figures7 and 14). At first this has 
the effect of reducing the lateral compositional non-uniformity. After 
45 seconds have elapsed since the termination of the impulse the 
composition nonuniformity is reduced to zero and subsequently 
increases in magnitude but has the opposite sense in comparison to 
the initial non-uniformity. This change in sense can be seen upon 
comparison of Figures 16 and 17. At approximately 260 seconds after 
the termination of the impulse the lateral segregation reaches a 
maximum value of 26%. It then decreases in value eventually reaching 
zero at 1017 seconds, changing sense and slowly increasing toward its 
initial steady level of 21.5% after more than 2000 seconds have 
elapsed. 

A shorter duration (10* 1 second) pulse resulted, 350 seconds 
after the termination of the pulse, in a maximum deviation of the 
lateral non-uniformity of only 5% from the initial steady level. 

The effects of two one second pulses separated by one second 
were also calculated. The magnitude of the pulses was 3x1 O' 3 g, and 
they were oriented parallel to the crystal interface. Their main effect 
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was to drive the lateral segregation from 21.5% (the initial value) to 
76% after 225 seconds. 

In addition double pulses were also examined. A pulse anti- 
parallel to the background steady acceleration followed by an equal but 
opposite pulse does not result in a "null" effect. While the flow 
generated by the first pulse is reversed by the second pulse there is a 
net flow following the termination of the second pulse. This flow is in 
the same sense as the initial steady flow and results in an increase of 
the lateral non-uniformity in composition to a maximum of 24% at 100 
seconds, whereupon it decays slowly to its initial value. 

4. Summary 

The salient results of our calculations for materials with 
properties and growth conditions similar to those listed In Table 2 can 
be summarized as follows: 

1) For a fixed growth rate, the amount of lateral non-uniformity in 
composition is very sensitive to the orientation of the steady 
component of the residual gravity vector. The worst case appears to 
be when the acceleration vector is parallel to the crystal interface. At 
growth rates on the order of microns per second, this orientation can 
lead to non-uniformities of 22% when the magnitude of the 
acceleration is V 2 ( 10)' 6 g. If, however the growth rate is lowered by 
an order of magnitude, the non-uniformity is reduced significantly 
(down to 4-5% in this case). 

2) A steady background level on the order of 10' 6 - 10* 5 g can be 
tolerated provided that the acceleration vector is aligned with the 
axis of the growth ampoule, and provided that no accelerations with 
frequencies less than 10‘ 2 Hz (and amplitudes of the order of the 
steady component) are present. 

3) The response of the solute field, and the lateral non-uniformity, 
to oscillatory accelerations varies from no response at all (at 
frequencies above 1 Hz with amplitudes below 10' 3 g) to a significant 
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response at 10-3 Hz at amplitudes on the order of 10' 6 g. In addition, 
additive effects were observed for combinations of a steady componen 
and low frequency components. These additive effects gave rise to 
significant lateral and longitudinal non-uniformities in concentration. 

4) The effects of impulse-type disturbances can be severe and can 
extend for a long time (on the order of 10» seconds) after the 
termination of the impulse. For example a pulse with a one secon 
du™uon or a combination of such pulses has a drastic effect on toe 
segregation levels at pulse amplitudes of 10-3 g. The nature of toe 
response depends on toe magnitude, direction and duration of toe 
impulse, and whether sequential opposing Impulses are involved A 
so-called ■ compensating" double pulse will not result *-o^tely 
offsetting effects. For the case we examined. however. toe resuWng 
compositional non-uniformity was not as severe as for sequent a 
pulses with the same orientation. Further investigation of more 
realistic impulses (g-JItter) Is necessary since toe response 
system appears to depend on toe nature of toe impulse and our results 
indicate that impulses appear to have important consequences 

transient behavior in crystal growth systems. 

It should be bome in mind that our calculations have only covered 
a small part of a large parameter space. In particular it should be 
noted that for a given level of residual acceleration the amount of 
lateral segregation can be expected to vary according to the 
magnitude of the Schmidt number and distribution coefficien . 
have also examined the effect of toe growth rate and have found that a 
reduction in growth rate by an order of magnitude wil result in a 
significant reduction in non-uniformity (for the ran S e of P ara ^ ete 
we have studied). This is consistent with the results of Adomato and 

Brown (1987). 
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Table l 

r ^lorotinn vector examined in this wor 



Steady 

^g 11^2(1 0)' 5 • 5^2(10)-6.’J2(10)- 6 . -J2(10)- 7 

S* 1 * g n cos(2no) n t): g(« - go * X gn cos(2*co„t, ; 

1- 10 -4. 10 -3. 10-2. 10-1. 1 10 HZ. 
uJonoH. S^2(10)-6. V2U01-6 


ft. 


impulse 

g(t) = gb t < 
g b = J2(10) 


ti, g(t) = gb + gl • ‘l < 1 < 
-6 ’ g , = 3 - J 2110)- 3 


t 2 . g(t) = gb ‘ > ‘2 



Table 2 

Thermo-physical properties characteristic of gallium-doped 
germanium ( after Chang & Brown 1983) and operating conditions 

Property 

Value 

Thermal conductivity of the melt 

0.17 WK’l cm'l 

Heat capacity of the melt 

0.39 Jg-lR-1 

Density of the melt 

5.6 g cm - 3 

Density of the solid 

5.6 g cm'3 

Kinematic viscosity of the melt (v) 

1.3 (10)-3 cm2 s -l 

Melting temperature (Tm) 

1231 K 

Solute diffusivity (D) 

1.3(10)-4 cm2 s -l 

Thermal diffusivity of the melt (k) 

1.3(10)‘l cm2s _ 1 

Segregation coefficient (k) 

0.1 

Thermal expansion coefficient (P) 

2.5 (10)' 4 K' 1 


Operating conditions 


Hot zone temperature (T H ) 

Distance between inlet and interface (L) 
Height of adiabatic zone 
Ampoule width (diameter) 

Translation (supply) rates (V M ) 

Associated dimensio nless parameters 

Prandtl number Pr = v/k 
Peclet number Pe = VmL/k 
S chmidt number Sc = v/D 
Density ratio a 


0.01 

5(10)"3 and 5(10)‘ 4 
10 
1.0 


1346 K 
1 cm 

2.5 mm 
1 cm 

6.5 pm s - 1,0.65 pm s‘l 


Compositional 


Table 3 

non-uniformity % [*] for computed 2-D (3D) 
steady state cases 


Residual 

Acceleration 

Magnitude 


io- 4 g 

10* 4 g 
V2(10' 5 ) g 
^2(10' 5 ) g 

10' 5 g 
10- 5 g 
5V2(10- 6 ) g 
V2(10" 6 ) g 
V2(10- 6 ) g 

■\/2(l0' 6 ) g 
V2(10- 6 ) g 
10' 6 g 
10' 7 g 


Orientation 

N c g 


\ 

I 

X 


4 


T 


I 

t 


l 


Ampoule Width 1 cm 2 cm 0.5 cm 
Growth Rate [pm s' 1 ] 

6.5 3.25 0.65 6.5 6.5 


36 

32 

110 

152 


(91) 

7.5 

57 

10 

22 

(26) 

4.0 

2.0 
0.7 
1.0 


4.6 0.7 


0.4 0.0 

0.5 0.2 
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3.8 


N _ normal vector to interface pointing 

e = unit vector parallel to g. w - non 

into melt. H = 1 cm for all cases 


figure captions ... . 

Figure 1 The temperature -composition phase diagram or a u 
binary system A-B. where B Is the dilute species. The 
compositions in the solid (S) and liquid (L) are 
respectively given by c s and c M . The equilibrium 
distribution coefficient, k, is given by the XY/XZ. 

Figure 2. Concentration profiles in the liquid and solid caused by 

segregation: 


Figure 3. 
Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


a) The zero growth rate (equilibrium) case. 

b) In this case the solute is rejected at finite 
velocities. This results in compositional gradients in both 

phases. 

The prototype directional solidification model. 

The dimensionless temperature field, 0. for all 2 D ca 
discussed in this paper. 

The steady flow field produced by a residual acceleration 
with a magnitude ^2(10)-5 g acting parallel to the crystal 
melt interface. The maximum speeds are approximate y 
twice the growth speed. 

The dimensionless solute field. C, associated with the flow 
depicted in Figure 5. For this case 152%. 

The steady flow field produced by a residual acceleration 
with a magnitude ^2(10)‘ 6 g acting parallel to the crystal 
melt interface. The maximum speeds are slightly greater 


than the growth speed. 

Figure 8. The dimensionless solute field. C. associated with the (low 
depicted in Figure 7. For this case $= 22%. 

Figure 9. The steady solute distribution (C s = c s /c~) in the crystal 
consequent to ^2 CIO) 5 g acceleration oriented parallel to 
the interface for the 3-D case. The cross section is taken 
perpendicular to the ampoule axis. \ - 91%. 


Figure 10. The steady solute distribution (Cs = Cs/Coo) in the crystal 
consequent to <2 (10) 6 g acceleration oriented parallel to 
the interface for the 3-D case. The cross section is taken 
perpendicular to the ampoule axis . £ = 26 %. 

Figure 11. Lateral non-uniformity in composition, plotted as a 
function of time for an oscillatory residual acceleration 
with a maximum magnitude of 3>/2(10)-6 g and a frequency 
of 10-3 Hz, acting parallel to the crystal-melt interface. 
The initial state was purely diffusive. 

Figure 12. Lateral non-uniformity in composition, plotted as a 
function of time for a residual acceleration consisting of a 
steady part with a magnitude of V2(10)* 6 g and an 
oscillatory part with a maximum magnitude of 3^2(10) 6 g 
and a frequency of 10*3 Hz, acting parallel to the crystal- 
melt interface. The calculation was started from a steady 
flow associated with a ^2(10)-6 g acceleration acting 
parallel to the interface. 

Figure 13 . Lateral non-uniformity in composition, plotted as a 
function of time for a multi-component disturbance 
consisting of a steady low g background plus two periodic 
components: g(t) = go + gi cos(27ilO- 3 t) + g2 cos(2jc10' t), 
where llg 0 M= ^2(10)* 6 , 1 1 gill = 3 V 2 ( 10)- 6 and 1 1 g 2 1 1 = 
3 ^ 2 ( 10 )" 3. The calculation was started from a steady flow 
associated with a V2(10)-6 g acceleration acting 
perpendicular to the interface. 

Figure 14 . The velocity field after a one second pulse of 3 ( 10)-3 g 
superimposed on a steady flow caused by a V 2 ( 10)' 6 g 
acceleration, both parallel to the crystal melt interface. 
The impulse, gi, is in the opposite direction to the 
background acceleration, gB- Note that the maximum 
velocity magnitudes are of the order 3 ( 10)' 2 cm s' 1 i.e. about 
500 times those in Figure 7 . 


Figure 15. 


Figure 16. 


Figure 17. 


Figure 18. 


Figure 19. 


Figure 20. 


Dimensionless solute field. C. immediately after the 

termination of the impulse. The compositional non- 
uniformity, is 21.5%. 

Dimensionless solute field. C. 31 seconds after the 

termination of the impulse. The compositional non- 
uniformity, is 6. 1%. 

Dimensionless solute field. C. 81 seconds after the 

termination of the impulse. The compositional non- 
uniformity. i is 11.4% and has changed sense. 

Dimensionless solute field, C. 431 seconds after the 
termination of the impulse. The compositional non- 
uniformity, Z>, is 21.7%. 

Dimensionless solute field. C. 881 seconds after the 
termination of the Impulse. The compositional non- 
uniformity. 5, is 4.5%. 

Dimensionless solute field. C. 1781 seconds after the 
termination of the impulse. The compositional non- 
uniformity, is 12.5%. 
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